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Abstract 

LISA (Laser Interferometer Space Antenna) is a joint mission of ESA and NASA which aims to be the first space-borne gravita- 
tional wave observatory. Due to the high complexity and technological challenges that LISA will face, ESA decided to launch a 
technological demonstrator, LISA Pathfinder. The payload of LISA Pathfinder is the so-called LISA Technology Package, and will 
be the highest sensitivity geodesic explorer flown to date. The LISA Technology Package is designed to measure relative acceler- 
ations between two test masses in nominal free fall (geodesic motion). The magnetic, thermal and radiation disturbances affecting 
the payload are monitored and dealt by the diagnostics subsystem. The diagnostics subsystem consists of several modules, and one 
of these is the magnetic diagnostics unit. Its main function is the assessment of differential acceleration noise between test masses 
D due to the magnetic effects. To do so, it has to determine the magnetic characteristics of the test masses, namely their magnetic 
remanences and susceptibilities. In this paper we show how this can be achieved to the desired accuracy. 
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1. Introduction 

LISA (Laser Interferometer Space Antenna) is a space mis- 
sion of NASA and ESA which aims at detecting low frequency 
gravitational waves. LISA will consist in a constellation of 
three spacecraft occupying the vertexes of an equilateral trian- 
gle of side 5 million kilometers. The barycenter of the con- 
stellation will orbit around the Sun following a quasi-circular 
orbit inclined 1° with respect to the ecliptic, and trailing the 
Earth by some 20°. Each of the spacecraft harbors two proof 
masses, carefully protected against external disturbances such 
as solar radiation pressure and charged particles, which ensures 
they are in nominal free-fall in the interplanetary gravitational 
field. Gravitational waves show up as differential accelerations 
between pairs of proof masses, and the main aim of LISA is 
to measure such acceleration using laser interferometry. The 
interested reader is referred to [ 1 1 and ^ for more extensive 
information. 

The technologies required for the LISA mission are many 
and challenging. This, coupled with the fact that some flight 
hardware cannot be tested on ground, led ESA to frame within 
its Scientific Program a technology demonstrator to test the re- 
quired critical technologies in a flight environment. Its launch 
is expected towards early 2014. The basic goal of LISA 
Pathfinder consists in measuring to the highest possible accu- 
racy the acceleration noise of two test masses 35 centimeters 
away, and to keep this noise below a certain limit OIU. The 
payload on board LISA Pathfinder is called the LISA Technol- 
ogy Package (LTP) fSjISl- Its main components are the two 
Gravitational Reference Sensors (GRSs) — shown as the two 
large vertical cylinders in figure [T] — and the Optical Metrol- 



ogy System (QMS) — which is an interferometer placed on the 
horizontal platform between them. The GRS consists in a set of 
electrodes aimed to determine the position of a test mass with 
respect to the spacecraft to nanometer precision, using capac- 
itance measurements. The QMS provides picometer precision 
measurements of the relative position of two test masses. Using 
these measurements, a set of micro-thrusters — Field Emission 
Electric Propulsion (FEEP) — relocate the spacecraft so that 
the test mass preserves a free fall motion. In the same figure [T[ 
the location of two induction coils can be observed. The coils 
will be used to generate a controlled magnetic field within the 
volume of the LTP Core Assembly (LCA) in such a way that the 
magnetic characteristics of the test masses can be determined. 
At the same time, the magnetic field inside the LCA will be 
measured by a set of four tri-axial fluxgate magnetometers (not 
shown in the figure). All these pieces of hardware constitute 
the Magnetic Diagnostic Subsystem. Finally, the Data Manage- 
ment Unit (DMU) is in charge of commanding and acquiring 
signals from all the mentioned subsystems. 

Magnetic noise in the LTP can be a significant part 
of the total readout noise, L2-10"''*ms"^Hz"'''^ out of 
S.OTO"''* m s"^ Hz"'^^ is the allocated maximum acceleration 
magnetic noise budget. This noise occurs because the resid- 
ual magnetization and susceptibility of the proof masses couple 
with a surrounding magnetic field, giving rise to a force Q: 
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and a torque: 

N = <MxB-Hrx [(M-V) B]) V 
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Figure 2: Coordinate reference frame of the experiment. 



Figure 1 : Ailistic view of the LISA Technology Package. The two towers are 
the Gravitational Reference Sensors. They are connected by the optical bench 
(grey), where the interferometer is located. The induction coils are located next 
to the two towers. 

In these expressions B is the magnetic field in the test mass, M 
stands for the magnetization of the test mass, V is the volume of 
the test mass, x is its magnetic susceptibility, /vq is the vacuum 
magnetic constant, and r denotes the distance to the center of 
the test mass. Finally, the notation (■ ■ ■) refers to the volume av- 
erage of the enclosed quantity. Thus, to estimate and ultimately 
subtract the acceleration noise due to the magnetic interactions, 
the magnetic properties of the test masses must be determined. 

In this paper we describe the experimental setup and the data 
analysis needed to infer the values of the magnetic properties of 
the test masses on board the LTP, and we assess the feasibility 
of obtaining the magnetic characteristics of the test masses with 
good accuracy. Specifically, we present a set of simulations 
aimed at evaluating the response of the LTP hardware (coils 
and test masses) and control architecture (drag free controllers 
and low frequency suspension controllers) when a controlled 
magnetic field is applied. We will show that using this pro- 
cedure, both the magnetization and the magnetic susceptibility 
of the proof masses can be determined to the desired accuracy. 
The paper is organized as follows. In section|2] we describe the 
main elements of the experimental setup. Section [3] is devoted 
to compute the forces and torques acting on the test masses. It 
follows section |4] where the different noise sources perturbing 
the experiment are presented. Finally, in section |5] we present 
our results, whereas in section |6] we summarize our main find- 
ings and we present our conclusions. 

2. Experiment description 

As mentioned, the basic approach to determine the magnetic 
properties of the test masses is to inject a controlled signal 
with the onboard coils and to study the dynamics of the proof 
masses. The two test masses are located at the center of each 
inertial sensor, and are the end mirrors of the OMS. In fact, 
one of the test masses will be the reference free floating body 
to perform the translation control of the spacecraft. The test 



masses are made of an alloy of Ft (27%) and Au (73%), their 
dimensions are 46 x 46 x 46 mm and their weight is 1 .95 kg. To 
comply with the top science requirements, the test masses must 
have certain properties. For the purpose of the present paper the 
two most important ones are the magnetic moment and the sus- 
ceptibility, which must be, respectively, |m| < 2.0 ■ 10"*^ A m^ 
and \x\ < 2.5 ■ 10"^ fSl. The volume of the test masses is 
V = 0.0046^ m^. The density of magnetic moment has to be 
then |M| < 9.451 ■ 10"^ A/m. 

The controlled magnetic field will be produced by the on- 
board coils, which are placed next to each of the GRS towers 
— see Fig. [T] The two circular induction coils are made of a 
titanium alloy (Tig AI4 V), have N - 2 400 windings of radius 
r - 56.5 mm. They are placed 85.5 mm away from the center 
of the respective test mass. The onboard coils are aligned with 
the X-axis of the test masses, thus, the magnetic field within the 
volume of the test masses has axial symmetry, see Fig. |2] If 
the current fed to the coils is /(f) = /q sin a>t, the resulting mag- 
netic field (and its gradient) will oscillate at the same frequency. 
Thus, when the coils are switched on the test masses rotate and 
are displaced from their equilibrium positions. Typical values 
of /q and (jj are, respectively, 1 mA and 1 mHz. 

Compared to other missions, LTP is a very flexible instru- 
ment in terms of the possible operation scenarios. Neverthe- 
less, we characterize the magnetic experiment for a fixed op- 
erating mode, the main science mode |9|. This mode is a full 
3-dimensional dynamical mode and it is schematically shown 
in Fig. |3] D is a dynamical matrix which represents the dy- 
namic response of the spacecraft and the test masses when they 
are aff'ected by specific forces (f). This block consists of an 
18 degree-of-freedom representation of the motion of these 3 
bodies. The diff'erential position of the test masses and their 
distance to the spacecraft are represented by x in this block- 
diagram. 

The LTP is endowed with two dififerent mechanisms to de- 
tect the motion and the actual position of the test masses. The 
first one is the interferometer, while the second is the electrode 
housing of each of the two gravitational reference sensors. The 
OMS in LTP is in charge of measuring the distance between one 
of the test masses and the optical bench, thus giving an abso- 
lute reference, and also the distance between both test masses. 
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Figure 3: Control system architecture of LISA Pathfinder. D stands for the 
dynamical matrix, S represents the sensing matrix of the inteiferometer, i.e. 
the matrix translating the position of the test mass, x, into the interferometer 
readout, o (On stands for the readout noise) . A represents the physics of the 
FEEPS and the electrostatic actuators, and finally C, is the controller matrix, 
implementing the drag free and low-frequency control loops. Oj represents the 
displacement guidance signals. a„ are the actuators noise and are the output 
forces of the actuators. f„,ag are the magnetic forces induced by the coils and f„ 
are the environment force noises disturbing the spacecraft. 



providing a differential reading. Due to its ability to perform 
wavefront sensing, the rotation angles of the test mass around 
the y- (r\) and z-axis ((f) can also be measured (the y- and z-axis 
in the test mass coordinate frame are the same axis represented 
in Fig. |2]but centered at the test mass). The displacements are 
expected to be measured with a picometer accuracy while the 
rotation angles can be measured with an accuracy of ~ 400 
nrad IfTOl [TTI . The electrode housing can also be used to deter- 
mine the position of the test masses. However, this mechanism 
only offers readings with nanometer precision. Consequently, 
for our application only the readings of the interferometer will 
be used. The physical model of these sensing mechanisms is 
included in S, the sensing matrix (figure [3]l. o„ is the readout 
noise of the different sensors and o the actual measure delivered 
to the controller 

These kinematic measurements are processed by the con- 
troller block (C) and a feedback action is produced on the LTP 
dynamics by the actuators (A): the micropropulsion thrusters 
and the electrostatic actuators. They produce an additional 
set of forces (fa) with the following objectives. The drag 
free controller acts on the spacecraft using the micropropulsion 
thrusters and forces it to follow test mass 1 . The electrostatic 
actuators act on test mass 2 using the low frequency suspension, 
a specific control loop of very low gain in the LTP measurement 
bandwidth, that allows to control in band the differential accel- 
eration between both test masses at the same time, avoiding 
secular drifts or stray motions of the second test mass. 

Finally, fn represents the force noise on the test masses, and 
fmag represents the forces acting on the test masses due to the 
magnetic field created by the coils. The complete control archi- 
tecture can be expressed by the following system of equations: 



O = D ' S f + On 

f = fmag + fn - A ■ C ■ (O H- Oi) - an 



(3) 



exception of Oi which represents the displacement guidance sig- 
nals of the experiment and an which are the actuators noise. 
Using Eqs. ([3]l, we can calculate the transfer function from the 
magnetic forces (f^ag) to the interferometer readings (o), which 
turns out to be: 



D 



fmag 1-HD-I-S-A-C 



(4) 



This transfer function characterizes the projection of the 
magnetic forces/torques into kinematic motion of the test 
masses. The controllers have been designed to deliver very sen- 
sitive readings of the differential motion of both test masses 
between 1 mHz and 30 mHz, the measurement bandwidth of 
the LTP mission. For simplicity, if we consider only the one- 
dimensional model: 
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where a>i and ai2 are the stiffness parameters coupling the mo- 
tion of each test mass to the motion of the spacecraft, dn is 
the interferometer channel crosscoupling. Cdf and Clfs are the 
drag free and the low frequency suspension controller transfer 
functions respectively and, finally, Afeep and Aea are the phys- 
ical models for the FEEP thrusters and the electrostatic actua- 
tors. 

3. Forces and torques 

Using Eq. ([T} and neglecting the environmental field (which 
is much smaller than the applied field, 500 nT), the x- 
component of the force acting on the test mass is 

F, = <M- VBo v> V sin Lut+— <Bo-VBo x) sin^ (6) 

/"() 



where Bq is the field produced by the coils. Thus, since 
sin^ cot = (1 - cos2(yf)/2, the linear acceleration of the test 
masses has two separate frequencies, one at oj and the other 
one at 2a), and also a DC component. The force on test mass 
1 is plotted in figure |4] The torque acting on the test mass also 
has a similar behavior. However, it must be noted that, because 
of the symmetry of the applied magnetic field, the torque only 
has one frequency component: 



N = <M X Bo + r X [(M-V) Bq]) V sin cjt 



(7) 



where all the symbols have been already defined with the only 



The resulting torques are displayed in figure|5] It is important 
to realize that only the y- and z-components of the torque can 
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Figure 4: The three components of the force on test mass 1 when coil 1 is on. 
These forces depend on the values of ^ and My, My and M-. For this specific 
example we have adopted M, = 16.4 ■ 10"' A/m, M,. = 9.1 ■ 10"' A/m, M- = 
-6.8 ■ 10"' A/m (a random orientation of the maximum M) and^ = 2.5 ■ 10"'. 
Iq is 1 mA and w is 1 mHz. The A-component of the force is shown as a solid 
black line, whereas the v-component is displayed as a dashed-dotted red line 
and the z-component is displayed as dashed green line. 



be measured with the interferometer, as A^v produces a rotation 
around the direction of the laser beam. Moreover, decomposing 
Eqs. (|6]l and (|7]i, it is easy to show that the x-component of the 
force on the test mass (which can be obtained from its displace- 
ment) and the y- and z-components of the torque (which can be 
obtained from the rotation angles) can be cast in the form: 



(8) 



where is a constant function, /v,^, /y,^ and oscillate 
at o) and fx2^ oscillates at Icj. Hence, A^y and Nr will be used 
to estimate and My, respectively, while the differential dis- 
placement of the test masses will be used to measure and;^'. 

4. Modelling of the noise sources 

The forces and torques shown in figures]?] and |5] correspond 
to an ideal case. However, the real forces and torques acting on 
the test masses will not be noise-free. Additionally, the outputs 
detected by the interferometer will also be aifected by several 
noise sources. Thus, to assess the feasibility of the experiment 
we need to model the noise sources. This is precisely the aim 
of this section. 



4.1. Magnetic hardware noise 

1 /2 

The stability of the magnetic field (S ^ ) produced by 
the coils at the position of the test masses and its gradient 
^^\1b idx^ must be, respectively, better than 5 nT Hz"'^^ and 
12 nT m"' Hz "'^^ [8| within the measurement bandwidth 
(1 mHz < / < 30 mHz). This can be translated into a require- 
ment on the stability of the injected current. It turns out that the 
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Figure 5: Torques on test mass 1 when coil 1 is on. The magnetic properties 
of the test mass are the same adopted in figure|4] The ,v-, y- and z-components 
of the torque are shown using black solid, red dashed-dotted and green dashed 
lines, respectively. 



requirement on the magnetic field gradient is the more demand- 
ing one, and using Ampere's law it is straightforward to show 

1/2 

that it is equivalent to a current fluctuation (S j ) requirement 
of 1 10 nA Hz"''^ within the measurement bandwidth. 
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Figure 6: Block diagram of the coil's current source. 

In figure |6] we show a block diagram of the different hard- 
ware stages used to produce the current that feeds the coils. 
This block diagram has two parts: a digital-to-analog process- 
ing stage and an analog processing stage. The first block con- 
tains two digital-to-analog converters (DAC), followed by tran- 
sresistance amplifiers (I/V) delivering low impedance output 
voltages. The first DAC sets the reference voltage of the sec- 
ond DAC and thus the peak amplitude of the sinusoidal current 
applied to the coil. The second one is configured with bipo- 
lar operation to generate the quantized signal with the previ- 
ously selected amplitude. In the analog processing block the 
signal is low-pass-filtered and amplified with a Howland cur- 
rent source |12|. Finally, a switch is used to select one of the 
three possible states: short-circuit, open-circuit or connected. 

1 /2 

The noise of the current source chain for a DC signal, S can 
be written as: 



S ^) - [G2GiS STAGE, + GiS STAGE2 + S STAGE,] '^^ (9) 

where 5 stage, is the noise density of the voltage reference 
and the first DAC, 5staget is the noise density of the tran- 
simpedance amplifier and the second DAC, S stage, is the noise 
density of the transimpedance amplifier, low-pass filter and 
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Howland current source, and Gi and G3 are, respectively, the 
gains in the the second and third stages. The noise density of 
the first stage can be obtained as a contribution of the first DAC, 
SoACi, and of the reference voltage, Sv,^!- However, the noise 
contribution of the first DAC is negligible with respect to that 
of the reference source, Vief . Thus, if the DAC works at its full 
scale, 5sTAGEi ~ ^v,^r Likewise, S'stagez can be computed as 
a combination of two different sources ^daCj and 5i/v, being 
SdaCj the noise density of the second DAC and S'l/v that of the 
transresistance amplifiers, and finally 5 stage, - 'S'lpf + 5'hcs, 
where 5lpf and 5hcs are, respectively, the noise densities of 
the low-pass filter and of the Howland current source. These 
noise sources have been measured for a DC current of 1 mA. 
The values obtained are 5stage, = 9.6nAHz"'^^, 5stage, = 
0.03 nAHz"'^^, and 5 stage, = 26.4 nAHz"'^^, respectively. 
The details of how these measurements were done are out of the 
scope of this paper and will be provided elsewhere. However, 
we mention that although the previous analysis has been per- 
formed for DC currents, when we operate at 1 mHz the domi- 
nant noise source is the quantization noise of the second DAC. 
This means that the total current noise can be modeled as: 



1/2 



(10) 



where the sole contributor is the quantization noise of the sec- 
ond stage, and thus the total noise assuming an uniform quan- 
tization and a signal amplitude greater than a quantization step 
is given by: 

2/0 1 



'^^ 2^^ ^/12T^ 



(11) 



where /s is the sampling frequency, A^b = 8 is the number of 
bits of the onboard hardware, and the rest of the symbols have 
been already defined. 

The DMU delivers the output at a rate of 1024 samples per 
cycle. Thus, the 1 mHz sinusoidal signal will be sampled at a 
frequency = 1.024 Hz. Since the highest sinusoidal current 
used in the experiment will be 1 mA, the highest noise will be 
2.22 //A Hz- This noise is 2 orders of magnitude larger than 
all other contributions and therefore will be the dominant noise. 
Figure|7]shows the current noise for coil 1 for a signal of 1 mHz 
and 1 mA (which are the nominal values of the experiment). As 
can be seen, the quantization noise level is above the DC current 
stability requirement for the mission (dashed green line) but, as 
it will be shown below, it is still sufficiently small to allow a 
reliable estimation of the magnetic properties of the test masses. 
It is also interesting to note that the quantization noise (thus the 
total noise) can be reduced when a 16-bit DAC is used — see 
Eq. ( [TT| . Finally, note that our theoretical results nicely match 
the experimental results. 

4.2. Other noise sources 

So far we have discussed the noise density of the onboard 
coils. However, the different electronic subsystems of the satel- 
lite produce magnetic fields, which are also noisy. The mag- 
netic noise produced by these subsystems (~ 50) has been mod- 
eled considering the fluctuating values of their magnetic mo- 
ments liT3L The expected displacement noise spectral density at 
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Figure 7: Noise density for a 1 mA sinusoidal current of 1 mHz nominal fre- 
quency. The dotted blue line shows the measured noise density, whereas the 
solid black line is our theoretical estimate (see text for details). The green 
dashed line represents the DC current fluctuations requirement. 



the positions of the test masses due to the magnetic field gener- 
ated by several electronic boxes is 3 • 10"" m Hz"'^^ at 1 mHz. 

Other noise sources that affect the measurements are the 
readout noise of the OMS, the noise induced by the GRS and 
that of the star tracker. As already mentioned, we only use the 
readings of the interferometer. However, the control chain uses 
the readings of the inertial sensors and of the star tracker to 
stabilize the drag free motion and the satellite attitude. There- 
fore, both the GRS and the star tracker affect the readings of 
the OMS, and must be considered. The dominant noise is the 
readout noise in displacement of the interferometer. We model 
it using a two pole/two zero noise shape filter Iil4ll : 



lA^OMsl = 



/ + 2- 10-V(27r) 
f + 2- 10-4/(27r) 



(12) 



At the positions of the test masses this results in a low frequency 
displacement noise of around 5 ■ 10 " m Hz '^^ at 1 mHz but 
it is the largest contributor at high frequency. 

The actuators are also a relevant noise source at low frequen- 
cies. This noise is due to the capacitive actuators and the FEEP 
thrusters which have to ensure the free fall motion of the test 
mass. The micro-propulsion system of LISA Pathfinder is com- 
posed of 12 FEEP thrusters. The force noise of each individual 
thruster is modeled as lfT4l : 



FEEP I 



/+10 
/+10- 



-2 



(13) 



Their contribution in the differential readout channel is around 
2 ■ 10 '° m Hz '^^ at 1 mHz, and it is somewhat larger on the 
absolute channel (which provides the distance from test mass 1 
to the spacecraft) where the FEEP thrusters are the main distur- 
bance source ID. 

Solar and infrared disturbances have been also modeled. So- 
lar disturbances are due to the solar flux impacting on the sur- 
faces of the spacecraft. Infrared disturbances are due to the 
infrared emission from the spacecraft external surfaces. These 
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Figure 8: Top panel: relative displacement of the test masses measured using 
the differential channel of the interferometer (solid black line). The dashed red 
line shows the absolute displacement with respect to the optical bench. Bottom 
panel: rotation about the y-axis — black solid line — and the z-axis — red line. 
Both panels show the response of test mass 1 when only coil 1 is fed with a 
current of 1 mA and 1 mHz. 

two noise sources are the most important contributors of dis- 
turbance on the spacecraft coordinates but they are highly at- 
tenuated by the control architecture and they turn out to be 
completely negligible in the interferometer readings. Finally, 
the test mass disturbance noise represent the internal distur- 
bances at the test mass coordinates. They are expected to be 
the most relevant sources at low frequency on the interferom- 
eter readings and their contribution has been estimated to be 
~ 3 ■ lO"'" m Hz"'^^ at 1 mHz. All these sources of noise affect 
the quality of the estimate of the magnetic properties of the test 
masses, as it will be shown below. 

5. Results 

We have estimated the magnetic properties of the test masses 
using the readings provided by the mission telemetry. The 
telemetry corresponding to the magnetic experiments will con- 
sist of the commands sent to the coils, the displacement read- 
ings of the interferometer (namely, the absolute and differential 
readings) and the wavefront rotation readings about the y- (77) 



and z- axes (0). For simplicity, we have assumed that the stiff- 
ness of the test masses, the actuators gains and the interferom- 
eter crosscoupling factors have already been determined f TSll . 

The simulated displacements and rotations measured by the 
onboard interferometer are displayed in figure |8] These dis- 
placements and rotations have been obtained integrating the 
equations of motion of a rigid solid, and including the drag 
free and low frequency controllers — see Sect. |2] Note that 
the problem has 18 degrees of freedom. In particular, each of 
the two test masses has 6 degrees of freedom, and the space- 
craft also has 6 degrees of freedom. The closed loop simu- 
lation is performed with appropriate simulation tools that will 
be used for mission operations [16. ,17J . As can be seen, the 
displacements of the test masses are below 8 nm, while in per- 
manent regime the corresponding rotations have amplitudes of 
~ 4 jurad. The very long transient of about 3 000 s of the rotation 
excusions — see the bottom panel of figure[8] — is due to the ef- 
fect of the low-frequency controller. This controller is designed 
to avoid drift excursions of the test masses with frequencies 
smaller than 1 mHz. Consequently, the transient is very long. 
A similar transient, although less evident, is present in the dif- 
ferential reading — see the top panel of figure [8] The reading 
of the displacement channel has two frequency components, w 
and 2(L), however, these components are difficult to see in the 
time series shown in figure [8] Additionally, these two compo- 
nents are not in phase with the forces shown in figure [4]because 
the LPF dynamics and the controllers introduce a phase delay 
to each of the two components. Finally, it is worth mention- 
ing that these displacements and rotations are within acceptable 
margins because they do not exceed the authority limits of the 
drag free and low-frequency controllers. 

To further illustrate the feasibility of the experiment, in fig- 
ure[9]we show the noise breakdown of the differential displace- 
ment reading of the interferometer. This figure has been ob- 
tained simulating the output of the entire instrument for each of 
the noise sources presented previously in Sect. [4] We simulated 
100000 seconds for each source of noise. Then, we performed 
the spectral estimation with a smoothed power spectral estima- 
tor based on the Welch estimator using a Blackman-Harris win- 
dow [1I8J . The time domain simulation and the spectral analysis 
have been performed using the LTPDA toolbox ifTTll . This is the 
data analysis tool that will be used for mission operations. As 
can be seen, in the frequency domain the signals at ll) and loj 
are clearly visible. Note that the most important contribution is 
that of the sensors noise, which is mainly characterized by the 
high-frequency noise of the interferometer, the contribution of 
the FEEP thrusters and the disturbance noise of the test masses. 
The environment noise, the magnetic hardware noise and of the 
solar and infrared emission contribution are totally negligible in 
the interferometer readings and do not represent any restriction 
in terms of parameter estimation quality. 

We have already shown that the displacements and rotations 
of the test masses can be detected even in the case in which all 
the noise sources are considered. Now the question to be an- 
swered is to which accuracy the magnetic properties of the test 
masses can be estimated. To this end we have used a classical 
linear least squares procedure |19|. The magnetic parameters 



6 



Table 1: Errors in the estimates of the magnetic properties. 





AM,. 


AM,. 


AM, 




No noise 




10-1^% 


10-1^% 


10-1^% 


Hardware noise 


0.13% 


0.08% 


0.09% 


0.12% 


Environmental noise 


0.12% 


0.26% 


0.24% 


0.10% 


Sensors noise 


0.87% 


0.97% 


1.05% 


1.01% 


Actuators noise 


0.96% 


0.99% 


1.25% 


1.17% 


Solar and infrared noise 


0.03% 


0.02% 


0.05% 


0.06% 


Test mass disturbance noise 


0.82% 


0.73% 


0.75% 


0.99% 


All sources 


1.15% 


1.53% 


1.72% 


1.25% 



are estimated in the following way. Let be the differential 
displacement signals from the interferometer, and Ry and R~ the 
rotation excursions around the y- and the z-axis, respectively, 
we write then: 



M, 
X 



R, 
R, 



My ■ r-^^ + rir. 



(14) 



where we have used Eq. (j8]l, d and r are the signals in displace- 
ment and rotation matched to the expected waveforms in u) and 
2(jL), as obtained from Eq. (jSjl, and n^, , and are the errors 
of the estimation model, namely, the displacement error, the 
T]- and the 0- error, respectively. Then the estimated magnetic 
properties of the test masses (M,, M,., M^ wdx) applying least 
square techniques are computed as: 



M, 
X 



■ rzj ■ Rz 



(15) 



As can be seen, the values of M, and;^' must be disentangled 
from D j because the dynamics of the test masses show two fre- 
quencies, while Ry and R, can be directly used to estimate the 
values of M, and M,.. 

We have examined the contribution of each of the noise 
sources in our estimation accuracy. The accuracy of the mea- 
surements of the magnetic properties of the test masses is 
mainly affected by the specific contribution of the noise source 
within the measurement bandwith. For instance, if some noise 
source has a relevant contribution around 1 mHz in the rotation 
signals, the estimating algorithm can not disentangle this con- 
tribution from that of the injected torque. In table [T| we list the 
accuracies of the estimated magnetic parameters obtained for 
each of the individual noise sources and that obtained when all 
the noise sources are present (last row). Rather naturally, the 




10 ■ 



10 



Frequency [Hz] 



10 



10 



Figure 9: Noise breakdown of each of the noise sources to the differential dis- 
placement signal received when a 1 mA and 1 mHz current is circulating in 
coil 1 . The black solid line shows the result of considering only the magnetic 
hardware noise. The red dashed line shows the result of considering only the 
environmental magnetic noise. The cyan dashed-dotted line only takes into ac- 
count the solar and infrared emission noise. The green dotted green, the sensors 
contribution. The solid blue line the actuators noise, whereas the black dashed 
line shows the result of considering the test mass noise. Finally, in magenta 
solid line the result of considering all the noise sources. 



results for and x shown in this table are closely related to 
the noise contributions shown in figure [9] — note that figure [9j 
only shows the noise breakdown for the differential displace- 
ment channel. The largest contribution to the error budget are, 
as expected, the actuators noise, the interferometer noise and 
the test mass disturbances. Nevertheless, the overall quality of 
the estimate is fairly good, 1.43% (mean square error of the rel- 
ative errors of all the estimated parameters). It is interesting to 
note as well that even if the rotations signals present signal to 
noise ratios around a factor of 3 smaller, we obtain errors of 
the same order of magnitude for the estimates of M,. and M,. 
This stems from the fact that the signals from which they are 
obtained must not be disentangled. 

6. Conclusions 

In this paper we have confirmed the feasibility of deriving the 
magnetic properties of the test masses of LISA Pathfinder. The 
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magnetic experiment is based on injecting controlled sinusoidal 
currents through the on-board coils and studying the dynamics 
of the test masses, as measured with the optical metrology sub- 
system. In our study we have performed numerical calculations 
that incorporate a full model of the dynamics of the test masses, 
realistic noise sources and up-to-date simulations of the inter- 
ferometer and inertial sensors. In particular, all the degrees of 
freedom of the test masses have been appropriately analyzed 
and we have fully taken into account the control architecture 
of LISA Pathfinder We have obtained that the displacements 
of the test masses along the x-direction are ~ 8 nm, while the 
rotation excursions are approximately 4/irad. These findings 
confirm that the magnetic experiment is within the authority 
margins of the drag free and low frequency suspension con- 
trollers. Consequently, any damage to the entire experiment 
when the coils are excited can be safely discarded. Moreover, 
we have shown as well that the displacement and rotation sig- 
nals can be processed and pipelined to an adequate estimation 
algorithm that allows to estimate both the magnetic moment and 
the magnetic susceptibility of the test masses to a good accu- 
racy. Specifically, assuming that the remnant magnetic moment 
is homogenous within the entire volume of the test masses, the 
estimates have errors below the 2% level. 
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